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ABSTRACT 


In  a  single-target  radar  environment,  matched  filters  provide  the  maximum  output 
signal-to-noise  ratio  for  target  detection  and  yield  the  minimum  mean-squared  er¬ 
ror  estimate  of  target  range.  In  a  multiple-target  environment,  the  sidelobes  of 
the  compressed  pulse  must  be  considered  in  the  system  design  because  of  the  like¬ 
lihood  of  false  alarms.  In  this  case,  the  signal  processor  uses  weighting  filters 
which  are  not  matched  to  the  transmitted  waveform.  In  this  report,  expressions 
for  the  mean-squared  range  estimation  error,  the  estimate  bias,  and  the  effects 
of  the  sidelobes  are  derived  in  terms  of  the  impulse  response  of  an  arbitrary  mis¬ 
matched  filter.  We  desire  to  find  that  impulse  response  which  leads  to  an  unbiased 
estimate  having  the  minimum  range  estimate  variance  subject  to  preassigned  reso¬ 
lution  {i.  e.,  sidelobe)  constraints.  This  optimization  problem  is  formulated  in 
state-space  in  which  the  optimal  control  law  is  sought.  Pontryagin’s  maximum 
principle  is  used  to  obtain  necessary  conditions  for  the  optimum  filter.  When  the 
sidelobe  constraints  are  neglected,  these  conditions  lead  to  the  matched  filter  solu¬ 
tion.  In  an  attempt  to  synthesize  the  optimal  filter  for  the  general  case,  we  set  up 
a  nonlinear  programming  problem  involving  the  set  of  unknown  Lagrange  multi¬ 
pliers.  This  should  be  a  computationally  easier  problem  to  solve  than  the  original 
variational  problem.  An  example  is  given  which  illustrates  the  methodology  for 
synthesizing  the  optimum  filter  when  the  class  of  admissible  controls  (i.  e.,  filters) 
is  restricted  by  physical  considerations.  It  is  in  this  case  that  the  real  power  of 
the  state-space  development  is  clearly  demonstrated. 
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OPTIMAL  MISMATCHED  FILTER  DESIGN 
FOR  RADAR  RANGING  AND  RESOLUTION 


I.  INTRODUCTION 

It  has  long  been  known  that  the  matched  filter  represents  the  "optimum"  means  of  processing 
data  to  obtain  estimates  of  target  range.  Optimality,  in  this  case,  means  that  the  estimates  have 
the  smallest  mean-squared  error  possible,  a  result  which  is  valid  only  when  the  signal-to-noise 
ratio  (SNR)  is  large.  In  many  applications,  target  resolvability  is  a  consideration  almost  as  im¬ 
portant  as  range  accuracy,  in  which  case  it  is  not  clear  that  the  matched  filter  is  the  best  re¬ 
ceiver  to  use.  In  fact,  for  some  waveforms  (for  example,  a  sinusoidal  pulse),  very  good  range 
accuracy  can  be  obtained  but  the  resolution  problem  is  significant  since  the  envelope  of  the 
matched  filter  output  signal  has  large  subsidiary  sidelobes. 

One  approach  to  this  problem  is  to  assume  that  the  receiver  is  a  matched  filter  and  then  try 

to  design  the  input  signal  which  will  produce  good  range  estimates  subject  to  constraints  on  the 

sidelobe  structure  of  the  compressed  pulse.  Algorithms  are  now  available  which  generate  the 
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solution  to  this  problem,  but  in  most  cases  the  resulting  waveform  is  quite  complicated,  making 
it  difficult  to  build  a  matched  filter. 

Another  approach  is  to  use  both  a  known  signal  which  can  be  transmitted  easily  and  a  mis- 
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matched  filter.  This  was  done  for  the  linear  FM  waveform  in  an  effort  to  reduce  the  sidelobes 
of  the  compressed  pulse,  but  at  the  expense  of  a  loss  in  range  accuracy.  Heretofore,  no  effort 
has  been  made  to  design  a  mismatched  filter  to  minimize  the  mean-squared  range  error  subject 
to  preassigned  constraints  on  the  sidelobe  structure. 

In  this  report,  we  assume  that  a  given  pulse  is  received  in  the  presence  of  additive  white 
Gaussian  noise  and  passed  through  a  filter  which  is  not  necessarily  matched  to  the  input  pulse. 

We  also  assume  that  the  range  estimate  is  made  by  locating  the  time  at  which  the  envelope  of 
the  filter  output  achieves  its  peak  value.  The  performance  of  this  estimation  scheme  has  been 
analyzed  previously  with  respect  to  measuring  the  loss  of  accuracy  and  detectability  due  to 
nonoptimum  filtering.  It  is  shown  that  the  inability  to  build  perfectly  matched  filters  does  re¬ 
sult  in  a  loss  of  accuracy.  However,  the  advantages  of  mismatching  with  respect  to  improving 
the  multiple-target  resolution  is  not  discussed.  In  general,  the  mismatched  filter  leads  to  bi¬ 
ased  estimates  of  the  target  range,  a  result  which  the  analysis  in  Ref.  3  fails  to  take  into  ac¬ 
count.  Since  the  estimate  bias  can  be  significant,  we  rederive  the  equations  describing  the 
performance  of  the  mismatched  filter  and  rectify  this  omission. 

In  addition,  the  resolution  properties  of  the  filter  are  derived  and  we  show  that  the  compos¬ 
ite  accuracy/resolution  performance  of  the  filter  depends  on  the  filter  impulse  response  and  its 
first  derivative.  We  then  formulate  an  optimum  control  problem  which  leads  to  the  filter  impulse 
response  resulting  in  the  best  range  accuracy  subject  to  preassigned  resolution  constraints. 
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An  attractive  feature  of  this  approach  to  the  design  problem  is  that  the  class  of  admissible 
impulse  responses  can  be  restricted  according  to  the  degree  of  complexity  one  is  willing  to  build 
into  the  receiver.  For  example,  if  tapped  delay  lines  are  to  be  used  in  the  realization,  then  the 
search  is  performed  over  the  tap  weights  and  spacings.  This  search  is  then  possible,  taking 
into  account  the  effects  of  tap  reflections  and  the  attenuation  characteristics  of  the  real  line. 

This  aspect  of  the  design  will  be  discussed  in  detail  in  a  subsequent  publication. 

Here,  we  concentrate  on  the  formulation  of  the  optimal  control  problem  which  leads  to  the 
best  mismatched  filter.  Using  the  maximum  principle,  we  derive  the  matched  filter  when  the 
sidelobes  are  ignored.  Then,  we  restrict  the  impulse  response  to  be  a  linear  combination  of 
known  basis  functions;  such  a  realization  is  useful  when  using  RC  lumped  parameter  networks. 
The  optimization  is  performed  over  the  weight  to  be  assigned  to  each  function.  In  the  Appendix, 
the  general  problem  with  sidelobes  is  analyzed  and,  using  the  maximum  principle,  we  are  able 
to  reduce  the  function  space  optimization  to  a  nonlinear  programming  problem  involving  the  un¬ 
known  multipliers. 

II.  SUBOPTIMAL  SIGNAL  PROCESSOR 

We  shall  assume  that  the  range  of  the  target  is  to  be  estimated  on  the  basis  of  the  received 
pulse 

r(t)  =  p(t  —  r  )  cos  (ojct  +  0)  +  n(t)  (II-l) 

where  tq  represents  the  true  time  delay,  0  is  an  unknown  phase  introduced  by  channel  disturb¬ 
ances,  and  n(t)  is  a  sample  function  of  a  zero  mean  Gaussian  random  process  with  covariance 
function  E[n(t)  n(t 1 )]  =  NQ6(t  —  t').  We  assume  that  this  signal  is  processed  by  the  receiver,  a 
block  diagram  of  which  is  shown  in  Fig.  1.  The  impulse  response  h(*)  is  arbitrary  except  for 
the  restriction  that  it  belongs  to  a  class  of  filters  H. 


2  tin  («ct) 


Fig.  1.  Signal  processor  for  range  estimation. 

Analyzing  the  receiver  operations  with  respect  to  the  waveform  in  Eq.  (II- 1),  we  see  that 


xQ(t)  =  [p(t  -  tq)  +  nc(t)]  cos  0  -  ng(t)  sin0 
xg(t)  =  [p(t  -  tq)  +  nc(t)]  sin0  +  ng(t)  cos  0 


(II-2a) 
( II— Zb) 
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where  nc(t)  and  ng(t)  are  the  quadrature  components  of  n(t)  which  can  be  expanded  as 


n(t)  =  n  (t)  cos  (a>  t  +  0)  —  n  (t)sin(u;  t  +  0) 
c  c  s  c 

where 

E[nc(t)  nc(tr)]  =  E[ng(t)  ng(t'))  =  NQ6(t  -t') 

E[nc(t)  ng(t')]  =  0  . 

The  filter  output  signals  are 

zc (t)  =  J  h(r  —  t)  xc(t)  dt 

z  (r)  =  C  Wt  -  t)  x  (t)  dt 

S  ^-00  S 

Substituting  Eqs.  ( II- 2)  into  (II-5),  we  obtain 

zq(t)  =  [y(T)  +  ^(t)]  cos0  —  V  3(r)  sin0 

zg(T)  =  [y(T)  +  t?c(t)]  sine  +  7?g(T)  cos  0 

where 


(II-3) 

(II-4a) 

(II-4b) 

(II-5a) 

(II  —  5b) 

(II-6a) 
( II- 6b) 


y(r )  =  r  h(r  -  t)  p(t  -  t  )  dt 

OO 


h(T  —  t)  nc(t)  dt 


(II-7a) 


(II- 7b) 


V  (t)  =  r  h(T  -  t)  n  (t)  dt 

S  oo  s 

After  summing  the  two  outputs  of  the  square-law  devices,  we  have  the  signal 
w(T)  =  (y(T)  +  7)c(r)]2  +  [r,g(T)]2  . 

The  receiver  declares  as  its  estimate  of  tq,  the  number  t  where 
oj(t)  =  max  oj(t) 

T 


( II— 7c) 


( II — 8 ) 


( II -9 ) 


In  order  for  the  proposed  estimation  scheme,  Eq.  (II-9),  to  lead  to  error-free  performance  in 
the  absence  of  noise,  it  is  necessary  that 


max  co( r )  =  c o(t  )  =  max  y  (t) 
t  °  r 


which  requires  that 

y*(r0)  =  o  . 

We  shall  see  that  this  condition  leads  to  unbiased  estimates  for  the  unknown  time  delay. 


(11-10) 


(11-11) 
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If  h(t)  =  p( —  t),  the  receiver  is  the  well-known  matched  filter  processor  and  its  perform- 
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ance  has  been  well-documented  in  the  literature.  The  performance  of  the  mismatched  filter, 
h(t)  7^  p ( — t ) #  for  estimating  target  range  was  given  some  attention  by  Hansen,  but  an  error  oc¬ 
curred  in  that  paper  because  it  was  assumed  that  the  processor  would  generate  an  unbiased  es¬ 
timate  for  any  filter  structure.  This  is  not  true,  in  general;  for  this  reason,  we  give  our  own 
derivation  of  the  range  accuracy  formula. 

In  addition,  we  would  like  to  point  out  the  fact  that  in  Ref.  3  the  motivation  was  to  determine 
the  loss  in  accuracy  as  a  result  of  using  an  imperfectly  constructed  matched  filter,  and  no  men¬ 
tion  was  made  of  the  idea  of  synthesizing  a  matched  filter  to  give  good  range  estimates  and  low 
sidelobes  simultaneously. 

To  proceed  with  the  analysis  of  the  mismatched  filter  (which  closely  follows  that  in  Ref.  4), 
we  assume  that  if  the  processor  is  to  be  any  good  then  the  estimate  r  should  be  close  to  t  ,  the 
true  value,  when  the  SNR  is  large.  Expanding  a ;(r)  about  tq,  we  obtain 


w(t)  =  u(to)  +  w'(t0)  (t  -  tq)  +  Y  w"(to)  (t  -  To)2  + 


(11-12) 


For  large  SNR,  the  higher  order  terms  can  be  neglected  and  the  maximum  of  cj(t)  occurs  at  the 
point  r  where  oj’It)  =  0.  Therefore, 


A 

T  —  T 


M'(To> 

W"(ro) 


(11-13) 


From  Eq.  (II-8),  we  see  that 

w'(t)  =  2[y(r)  +  iJc(r)]  [y'(r)  +  tj^(t)]  +  2 rj g(r)  tj^(t) 


(11-14) 


w"(t)  =  2[y'(r)  +  t?c(t))2  +  2[y(T)  +  ry c(t)]  [y"(r)  +  7?y(T)] 


+  2[r)^(r)]2  +  2r,s(r)  v's'(r)  .  (11-15) 

To  find  an  expression  of  (r  —  tq)  which  is  first  order  in  the  noise,  we  express  Eq.  (11-14)  to  first 
order  so  that 

{«'(Tol=y(ro)y'(rol+y(T0|il'(To)+y'(To)iic(r)  .  (11-16) 

We  shall  assume  that  the  filter  has  been  designed  for  perfect  noise-free  performance,  so  that 
y'(To)  =  0,  which  follows  from  our  earlier  discussion.  Then,  Eq.  (11-16)  becomes 


}u'(T0)=y<Tc)^<ro)  .  (n-17) 

Since  this  is  first  order  in  the  noise,  it  suffices  to  express  wm(to)  to  zero  order  in  the  noise. 
From  Eq.  (11-15),  we  find  that 

J  m"(t0)  =  y(T0)  y"(ro)  (11-18) 


where  we  have  once  again  made  use  of  the  fact  that  y'(ro)  =  0-  Then,  the  error  in  the  range  es¬ 
timate  is  conveniently  expressed  as 


A 

T 


To  -y ”(t  ) 


rj'  (r  ) 
'c  o7 


(11-19) 
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Referring  to  Eq.  (II-7b),  we  see  that 


J’c(To)=Lh'(To-t)  ”c(t)  dt 


which  is  a  Gaussian  random  variable  with  mean  zero  and  variance 


<JZ  =  N  f  [h’(-t)]2dt 

°  V_  OO 


From  Eq.  (II-7a),  we  see  that 


Joo  s\  OO 

h(T  -  t)  p(t  -  t  )  dt  =  \  h( — t )  p(t  +  t  -  t  )  dt 

.00  °  ^Loo  ° 


(11-20) 


(11-21) 


(11-22) 


and,  therefore. 


/V» 

y ' <T0)  =  \  h(-t)  p'(t)  dt  (II-23a) 

V—  oc 


y"(T  )  =  f  h( —  t)  p"(t)  dt  .  ( II—  2 3b) 

U  OO 


Relating  Eqs.  (11-20)  and  (11-23)  to  (11-19),  we  conclude  that  for  large  SNR  the  proposed  estima¬ 
tion  scheme  produces  range  estimates  which  are  unbiased  estimates  of  the  true  parameter  value 

t  ,  and  have  a  Gaussian  distribution  about  t  .  The  mean-squared  estimation  error  is 
o  o 


E(t  —  t  )2  =  N 
o  o 


/.”  th'(-t)]2  dt 
t /_°L  h(-t)  p"(t)  dt]2 


(11-24) 


where  E  denotes  the  ensemble  average.  This  result  is  valid  for  the  class  of  filters  for  which 


C  h( —  t)  p'(t)  dt  =  0  .  (11-25) 


In  fact,  the  proposed  processor  is  meaningful  only  for  this  class  of  filters. 

If  the  filter  is  matched  to  the  signal,  h(t)  =  p( —  t)  and,  in  this  case,  Eq.  (11-25)  is  always 
satisfied.  It  is  then  easy  to  show  that  Eq.  (11-24)  becomes 


E(?-to)2  = 


N 


/_°L  [p'<t)r  dt 


(11-26) 


which  is  the  classical  result,  and  provides  a  check  on  our  analysis  to  this  point. 

Thus  far,  the  analysis  has  been  restricted  to  the  large  SNR  range  accuracy  performance 

of  the  receiver.  The  multiple-target  resolution  capabilities  can  be  taken  into  account  by  first 
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observing  that  the  envelope  y  (t)  will  have,  in  addition  to  a  maximum  at  tq,  subsidiary  peaks  — 
called  sidelobes  at  other  r  values.  If  these  sidelobes  are  the  same  order  of  magnitude  as 
y2(r  ),  it  will  not  be  possible  to  distinguish  between  two  distinct  but  neighboring  targets.  There¬ 
fore,  it  is  desirable  to  design  the  filter  to  make  these  sidelobes  small  with  respect  to  the  mag¬ 
nitude  of  the  central  lobe  at  r  .  This  can  be  done  by  requiring  that  the  constraint 

y2(*r)4  e(T)  y 2(tq) 
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be  satisfied,  where  €(r)  represents  the  sidelobe  constraint  function  and  is  chosen  to  combat  the 
particular  clutter  distribution  under  consideration.  Since 


y(r) 


h(r  —  t)  p(t  —  r  )  dt 
h(-t)  p(t  +  r  —  tq)  dt 


(11-27) 


then  y(TQ)  is  independent  of  and  the  sidelobes  can  be  kept  low  by  requiring  that 


If. h'-’ 


II 


t)  p(t  +  r)  dt  -  e(r)  h( —  t )  p(t)  dt  4  0 


(11-28) 


In  practice,  h(-)  and  p( • )  are  band-limited  functions,  in  some  sense,  and  the  continuum  of  con¬ 
straints  can  be  replaced  by  the  finite  number  of  constraints 


in 


h( — t )  p(t  +  7\)  dt 


-  h( — t)  p(t)  dtj  4  0  j  =  1,  2,  .  . .  ,  n 


(II- 29) 


The  foregoing  analysis  leads  to  the  following  filter  design  problem:  From  the  class  admis¬ 
sible  filters  H,  we  wish  to  find  that  filter  which  minimizes  the  quantity 


[h'<-t)]Z  dt 


(II-30a) 


[ h(-t)  p"(t)  dt]2 
subject  to  the  "zero-bias”  constraint 

f  h( —  t)  p'(t)  dt  =  0  (II-30b) 

J-  oo 

and  multiple-target  resolution  constraints 

[j  h(-t)  p(t  +  r.)  dtj 2  -e.  [j  h(— t)  p(t)  dtj  40  j  =  1,  2 . n  .  (II-30c) 


In  Sec.  Ill,  we  use  state-space  techniques  to  formulate  this  design  problem  as  an  optimal 
control  problem  in  state-space.  Then,  by  using  the  maximum  principle,  we  can  derive  the  con¬ 
ditions  necessary  for  optimality. 


III.  STATE-SPACE  FORMULATION  OF  THE  DESIGN  PROBLEM 

Here,  we  shall  formulate  the  mismatched  filter  design  problem  using  state-variable  tech¬ 
niques  so  that  the  theory  of  optimal  control  can  be  used  to  synthesize  the  optimum  filter.  In  or¬ 
der  to  do  this,  we  first  assume  that  the  filter  impulse  response  is  of  finite  duration,  T  seconds, 
and  is  identically  zero  for  t  i  [0,  T].  In  addition,  we  assume  that  the  transmitted  pulse  p(t)  is 
also  time  limited  to  the  interval  [0,  T),  an  assumption  which  will  always  be  satisfied  in  practice. 
Finally,  we  point  out  the  fact  that  in  the  analysis  of  Sec.  II  the  performance  depended  on  that 
part  of  the  impulse  response  for  t  <  0.  In  order  to  guarantee  the  realizability  of  the  optimum 
filter,  we  need  only  replace  h( — t)  in  all  our  equations  by  h(T  —  t). 
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We  define  the  control  function  u(t)  as  the  first  derivative  of  the  impulse  response,  while  the 
impulse  response  itself  is  set  equal  to  the  first  component  of  the  state  vector;  i.  e., 

u(t)  =  h'(T  —  t)  (III-l) 

x1(t)=h(T-t)  (III-2) 


which  leads  to  the  state  equation 

x^tjr-utt)  { III — 3 ) 

and,  since  the  filter  is  initially  at  rest  when  the  pulse  arrives,  we  set 

Xj(0)  =  0  .  (III-4) 

In  addition,  we  define  four  more  state  variables  according  to  the  state  equations 

x2(t)  =  h(T  -t)  p(t)  ( III- 5a) 

x3(t)  =  h(T  -  t)  p'(t)  (III- 5b) 

x4(t)  =  h(T -t)  p"(t)  ( III-5c) 

x5(t)  =  [h'(T  -t))2  ( III— 5d ) 

each  having  the  initial  conditions 

x2(0)  =  x3(0)  =  x4(0)  =  x5(0)  =  0  .  ( III- 6 ) 


It  is  obvious  that 


x2m  = 

fT  h(T  -t)  p(t)  dt 

J0 

(III- 7a) 

X3(T)  - , 

f  h(T  -  t)  p'(t)  dt 

Jo 

(III- 7b) 

x4m  = 

h(T  -  t)  p"(t)  dt 

Jo 

( III— 7c) 

x5m  = 

fT  [h'(T  -  t))2  dt  . 

Jo 

( III- 7d ) 

Furthermore,  by  substituting  Eqs.  ( III- 1 )  and  (III-2)  into  (III-5)  we  obtain  the  following  set  of 
first-order  differential  equations 


xd(t)  =  — u(t) 

x1(0)  =  0 

(III-8a) 

x2(t)  =  p(t)  x^t) 

x2(0)  =  0 

( III— 8b) 

x3(t)  =  p'(t)  x1(t) 

x3(0)  =  0 

( III— 8c ) 

x4(t)  =  p"(t)  x^t) 

x4(0)  =  0 

(III— 8d) 

i5(t)  =u2(t) 

x5(0)  =  0  . 

( III- 8e) 
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The  zero  bias  constraint,  Eq.  ( II- 30b),  requires  that 


x3(T)  =  0 


(III-9) 


which  follows  from  the  definition  Eqs.  (Ill- 5b)  and  ( III- 7b) .  Subject  to  this  constraint,  the  best 
range  estimate  is  obtained  by  minimizing  Eq.  (II- 30a),  which  is  equivalent  to  minimizing 


-x^(T)/x5(T) 


(III-10) 


which  also  follows  from  the  above  definitions. 

Finally,  the  multiple-target  resolution  requirements  are  accounted  for  by  defining  the  state 
variables 


y.(t)  =  h(T-t)  p(t  +  r.)  j  =  1,2, 

with  yj(0)  =  0  for  all  j-  Then  it  is  clear  that 


,n 


y.(T)  =  fT  h(T  -  t)  p(t  +  t  )  dt 
J  JQ  J 


and  for  good  resolution,  we  therefore  require  that 

j  =  1,  2,  .  . 


y.2(T)-e.x2(T)4  0 


(III-ll) 


( III— 12) 


( III— 13) 


which  follows  from  Eq.  ( II—  3 0c )  and  the  above  definitions.  Therefore,  the  filter  design  problem 
is  equivalent  to  the  following  optimal  control  problem.  Find  the  control  function  u(t)  and  state 
vector  [X| ( t), ....  x5(t),  y1(t), ....  yn(t)J  which 

(a)  satisfy  the  differential  equations 


xd(t) 

=  -u(t) 

x2(t) 

=  p(t)  x1  ( t) 

x3(t) 

=  p'(t)  x1(t) 

x4(t) 

=  p"(t)  x^t) 

x5(t) 

=  u2(t) 

^(t) 

=  p(t  +  Tj)  X1(t) 

j  =  1.2, 

(b)  satisfy  the  initial  conditions 

x1(0)  =  ...  =  x5(0)  =  0 

yt(0)  =  . . .  =  yn(0)  =  o 

(c)  satisfy  the  terminal  conditions 

x3(T)  =  0 

y2(T)- e  x|(T)^  0  j  =  1,2,- 
J  J 
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(d)  minimize  the  functional 
-x42(T)/x5(T)  . 

The  optimum  control  function  is  to  be  selected  from  some  class  of  admissible  controls  U. 
This  class  is  directly  related  to  that  of  admissible  impulse  responses  H.  So  far,  we  have  im¬ 
plicitly  assumed  that  if  h  e  H,  it  is  at  least  once  differentiable  and,  therefore,  u  c  U  must  be 
piecewise  continuous.  Additional  restrictions  on  U  can  be  imposed  by  taking  into  consideration 
the  specific  application  to  which  the  filter  is  to  be  used.  For  example,  in  a  radar  application, 
it  might  be  necessary  to  build  the  filter  using  a  tapped  delay  line.  The  optimization  would  then 
be  done  on  the  set  of  tap  weights,  which  merely  requires  imbedding  the  structure  of  the  delay 
line  into  the  above  formulation.  Even  the  physical  properties  of  the  line,  such  as  tap  reflections 
and  line  attenuation,  can  be  handled  in  this  manner.  This  represents  the  real  power  of  the  state- 
space  approach,  since  no  matter  how  complicated  the  system  may  become,  efficient  algorithms 
can  be  brought  to  bear  on  the  problem.  We  are  currently  investigating  this  aspect  of  the  design 
for  linear  FM  pulses,  and  will  report  our  findings  thoroughly  in  a  future  publication. 

IV.  EXAMPLES  OF  CONDITIONS  NECESSARY  FOR  OPTIMALITY 

We  shall  now  illustrate  the  control  theoretical  techniques  which  must  be  used  to  solve  for 
the  optimum  control  analytically.  In  the  first  example,  we  shall  neglect  the  resolution  con¬ 
straints  and  allow  the  class  U  to  be  all  piecewise  continuous  functions.  Of  course,  the  optimum 
filter  will  be  the  matched  filter.  The  result  is  interesting,  however,  since  it  illustrates  the 
fact  that  variational  methods  applied  to  the  estimate  variance  for  mismatched  filters  lead  to 
the  matched  filter.  Heretofore,  this  result  has  been  obtained  only  by  using  statistical  methods 
as  in  Ref.  4.  Control  techniques  have  been  used  previously  to  derive  the  matched  filter,5  but  in 
those  papers,  usually  only  an  artificially  defined  SNR  is  optimized.  We  have  therefore  obtained 
a  rather  ideal  blend  of  the  control  and  communication  theoretical  techniques. 

In  the  second  example,  we  restrict  the  class  of  filters  H  to  be  linear  combinations  of  orthog¬ 
onal  basis  functions.  The  optimum  weighting  for  each  component  is  obtained  analytically,  again 
assuming  no  sidelobe  constraints.  This  result  is  useful  in  that  it  shows  how  the  method  applies 
to  a  smaller  class  of  filters  and,  in  addition,  functions  can  be  chosen  which  lead  to  a  convenient 
RC  realization  of  the  optimum  filter.  Although  such  a  filter  cannot  perform  as  well  as  a  matched 
filter,  it  represents  the  best  possible  RC  approximation  to  the  matched  filter,  in  the  sense  of 
minimum  estimate  variance.  This,  then,  is  certainly  a  better  approach  than  using  cut-and-try 
to  get  close  to  matched  filter  performance  as  was  done  essentially  in  Ref.  3. 

The  third  example  is  a  repeat  of  the  first  except  that  the  sidelobe  constraints  are  considered. 
Since  most  of  the  manipulations  are  similar  to  those  used  in  the  first  two  examples  above,  the 
details  are  presented  in  the  Appendix.  In  this  case,  it  is  not  possible  to  solve  for  the  optimum 
control  analytically.  However,  we  show  that  the  variational  problem  can  be  reduced  to  a  non¬ 
linear  programming  problem  involving  the  unknown  multipliers.  This  should  be  a  simpler  prob¬ 
lem  to  solve  numerically  than  the  first.  We  now  consider  each  of  the  examples  in  detail. 

Example  1:  No  resolution  constraints,  arbitrary  U 

In  this  case  the  state  equations  are  simply 

3^(0= -u(t)  xd(0)  =  0  (IV-la) 
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x2(t)  =  p(t)  x^t)  x2(0)  =  0  (IV-lb) 

x3(t)  =  p'(t)  xd(t)  x3(0)  =  0  (IV-lc) 

x4(t)  =  p"(t)  x3(t)  x4(0)  =  0  (IV- Id) 

x5(t)  =  u2(t)  x5(0)  =  0  .  (IV-le) 

Defining  the  functions 

eQ[x(T)]  =  -x42(T)/x5(T)  (IV-2a) 

ejxfT)]  =  x3(T)  (IV-2b) 


we  want  to  minimize  0Q  subject  to  the  constraint  9^=0.  The  Hamiltonian  for  this  problem  is 
H(x,  X,  u)  =  — X^(t)  u(t)  +  X 2 (t)  p(t)  x3(t)  +  X 3<t)  p'(t)  X| (t) 


+  X4(t)  p"(t)  x4(t)  +  X5(t)  u2(t)  (IV-3) 

where  the  costate  variables  Mt)  satisfy  the  equations 

x  .(t)  =  -P-  j  =  1 . 5  .  (IV  -4) 

3  j 

Therefore, 

X3(t)  =  -p(t)  X2(t)  -  p'(t)  X3(t)  -  p"(t)  X4(t)  (IV-5a) 

X2(t)  =  0  (IV-5b) 

X3(t)  =  0  (IV-5c) 

X4(t)  =  0  (IV-5d) 

X5(t)  =  0  .  (IV-5e) 


The  terminal  values  of  the  costate  variables  are  obtained  from 

i  a©:[x(T)] 

YTI  ■  2  “i  ”axTfr 

i=0  3 

and,  therefore. 


Ad(T) 

=  0 

(IV- 7a) 

x2(T) 

=  0 

(IV- 7b) 

x3(T) 

=  or4 

(IV-7c) 

x4(T) 

=  -2aQx4(T)/x5(T) 

(IV-7d) 

X5(T) 

=  aox2(T)/x52(T)  . 

(IV-7e) 
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Combining  Eqs.(IV-5)  and  (IV-7),  we  get 


*2(t)  s  0  (IV-8a) 

X3(t)  =  »1  =  X3  (IV-8b) 

X4(t)  =  -2aox4(T)/x5(T)  =  X4  (IV-8c) 

X5(t)  =  qox42(T)/x52(T)  ^X5  (IV— 8d ) 

and  the  equation  for  X ^ ( t)  is  simply 

X^t)  =  -X3p,(t)  -  X4pM(t)  (IV -9) 

with  X |(T)  =  0.  This  equation  is  easily  integrated  to  give 

X^t)  =  X3[p(T)  -  p(t)]  +  X4[p'(T)  -  p'(t)]  .  (IV-10) 

We  shall  assume  that 

p(  0)  =  p(T)  =  0  (IV-lla) 

p*  ( 0 )  =  p'(T)  =  0  (IV-llb) 

so  that 

X^t)  =  -X3p(t)  -X4p'(t)  .  (IV- 12) 

The  maximum  principle^  states  that,  for  the  optimum  control,  X(t)  ^  0  and  0,  and  that  the 

Hamiltonian  is  maximized  by  the  optimal  control.  This  is  equivalent  to  maximizing 

h(u)  =  — X^(t)  u(t)  +  X^u^(t)  .  ( IV  —13) 


The  maximization  is  performed  over  the  set  of  u  e  U.  At  this  point,  the  properties  of  the  real 
filter  are  taken  into  account  by  restrictions  imposed  on  U.  For  the  first  example,  U  is  arbi¬ 
trary  and  the  methods  of  ordinary  calculus  can  be  used  in  the  maximization.  There  are  two 
cases  to  consider  at  this  point:  otQ  =  0,  and  <  0.  If  =  0,  then  from  Eq.  (IV-8),  X 4( t )  =  0, 

X^(t)  =  0,  and  in  Eq.  ( IV -12),  X ^ (t )  reduces  to 

X  !  (t )  =  —  X  3p(  t ) 

If  X3  =  0  or  X ^ ( t )  =  0  for  some  t,  then  X(t)  =  0  and  this  contradicts  the  necessary  condition  for 
optimality.  Therefore,  X3  ^  0  and  maximizing  Eq.  ( IV- 13)  means  maximizing 

h(u)  =  X 3p(t)  u(t) 


which  is  meaningless.  We  therefore  conclude  that  o tQ  <  0,  in  which  case  X,.  <  0  from  Eq.  (IV-8d). 
Therefore,  Eq.  ( IV -13)  has  the  well-defined  maximum 

u(t)  =  X1(t)/2X5  . 

Using  Eq.  (IV- 12),  we  conclude  that  the  optimum  control  is  of  the  form 


“(t>  =-2xT  P(t)“ 


2x7  p'<« 


(IV-14) 
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Thus,  the  optimum  filter  is  of  the  form 


X3  P*  X4 

l(t)  =  JQ  P(ff)  dCT  +  2X7  P(t)  • 


(IV-15) 


The  zero  bias  constraint  requires  that  x^(T)  =  0,  which,  from  Eq.  (IV-lc),  is  equivalent  to 


rT 

\  p'(t)  x.(t)  dt  =  0 

1 


(IV-16) 


Applied  to  Eq.  (IV-15),  we  get 


^  p* (t)  p(cr)  daj  dt  +  J  p* (t)  p(t)  dt  =  0 

Using  integration  by  parts  and  the  conditions  in  Eqs.  (IV-11),  we  get 


_J!  f' 

J, 


p  (t)  dt  =  0 


5  ^o 

Since  the  integral  represents  the  energy  in  the  transmitted  signal,  which  is  nonzero,  then  it  is 
necessary  that  X^  =  0.  By  using  this  result  in  Eqs.  (IV-14)  and  (IV-15),  the  optimum  control  and 
filter  are 


u(t)  =~2X~  P,(t> 


(IV-17a) 


xl(t*  =  2X  p(t) 

5 

The  minimum  cost  is  — x^(T)/x^(T)  where 


x4(T)  =  J  p"(t)  x1(t)  dt 

x4  rT 

=  2X74  p"(t)p(t) 


dt 


Xg(T)  =  JT  u2(t)  dt 


(IV-17b) 


(IV-18a) 


X4  pT 

-TJ  [P' 


(t)]2  dt 


(IV-18b) 


Using  integration  by  parts  and  Eqs.  (IV-11),  we  can  show  that 
r>T  r*T 


^  p”(t)  p(t)  dt  =  -  j*  [p1  (t)]2  dt 


and  therefore  the  minimum  cost  is 

4> 


-  JT  ip'(u]! 


dt 


(IV-19) 


(IV-20) 
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Since  this  is  independent  of  the  unknown  multiplier  X^/ZX^,  we  might  as  well  set  it  equal  to  unity, 
in  which  case 


u(t)  =  -p'(t) 


(IV-21a) 

(IV-21b) 


x\(t)  =  p(t)  . 

However,  we  defined  ( t)  =  h(T  —  t).  Therefore,  the  optimum  impulse  response  is 

h(t)  =  p(T  -  t)  ( IV-22) 

which  defines  the  matched  filter. 

Example  2:  No  resolution  constraints,  restricted  U 

The  only  difference  between  this  and  the  preceding  case  is  the  class  of  functions  U  over 
which  the  Hamiltonian  is  to  be  maximized.  Recall  from  Eq.  (IV - 1 3)  that  the  optimum  control  in 
the  set  U  maximized  the  function 


h[u(t),  t]  =  —X  1<t)  u(t)  +  X5u  (t) 

where 

xd(t)  =  — X3p(t)  -X4p'(t)  . 

In  this  case,  we  set 


U  = 


u:u(t)  =  ]] 

j  =  l 


( IV-23) 


(IV-24) 


( 1V-25) 


where  are  orthonormal,  i.  e., 

J  *.(t>  *k(t>  * « «.k 


for  all  j,  k 


( IV-26) 


We  shall  use  vector  notation  with  (•)t  denoting  the  vector  transpose.  Then,  if  u  e  U  maximizes 
Eq.  (IV-23)  for  each  t  (i.  e.,  u  is  the  optimal  control),  then 

h[u(t),t]4  h[u(t),t]  for  all  t  e  [0,  T]  (IV-27) 

and,  therefore, 

PT  PT  /v 

\  h[u(t),  t]  dt<$  \  h[u(t),t]dt  .  ( IV-28) 

Jo  Jo 

Since  u  e  U,  then  u(t)  =  ^(t)  c.  Using  this  fact  and  Eq.  (IV-23)  in  Eq.  (IV-28),  we  find  that 
[-  JT  Xd(t)  /(t)  dtj  c  +  X5ct  |  J  0(t)  ^(t)  dt  j  c 


|-  J  X1(t)  ^(t)  dt  J  £  +  X^c"  |  J  0(t)  ^(t)  dtj  c 


(IV-29) 
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where  c  leads  to  the  optimal  control  u.  Since  the  basis  functions  were  orthonormal,  the  matrix 


4  ( t)  4>t 


(t)  dt  =  I 


C¥  =  -  J  X1(t)  lp 


tf>(t)  dt 


which  is  the  m  X  m  identity  matrix.  Letting 

"»T 

Jo 

we  conclude  that  the  optimum  set  of  coefficients  c  must  maximize  the  function 
lie)  =  a}c  +  X^c 

We  have  already  shown  that  X^  <  0;  therefore,  the  c  which  maximizes  Eq.  ( IV -31)  is 
c=-a/Z\s  . 


But,  from  Eqs.  (IV-30)  and  (IV-24), 


A  rT  A  pT 

£=2^jo  p(t)  ^(t)  dt  +  ia7  i  p,(t)  a 


\  fT 

5  4) 


ip(t)  dt 


We  let 


p'(t)  t)  dt 


(IV-30) 


(IV-31) 


(IV-32) 


(IV-33) 


(IV-34a) 


b  =  ^  p(t)0(t)dt  (IV-34b) 

*~o 

and,  therefore, 

•  HV-351 

Substituting  this  back  into  the  defining  equation  for  u  e  U,  we  see  that  the  optimal  control  is 

(IV-36) 

The  corresponding  state  variable  x^is 

\A  X_ 

x^t)  =  “  £(t)  2X  ^(t)  ^  (IV-37) 

5  5 

where 

nt  t 

t(a)  dcr  .  (IV-38) 


u(t)  =  ^(t)  a  +  ^(t)  b  . 

5  5 


fl(t)  =  J  4>}( 


It  is  convenient  to  define  new  functions 
Xa<t)  =  ^(t)  a 

Xb<t)  =  f^t)  b 


(IV-39a) 

(IV-39b) 
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(IV-40a) 


$a(t)  =  ^(t)  a 


ii 

(IV-40b) 

and  variables 

a1  =  A4/2X5 

(IV-41a) 

a2  =  A3/2X5  . 

(IV-41b) 

Notice  that  while  a ^  and  a ^  are  as  yet  unknown,  the  functions  xa. 

Xb.  *a.  and  4>b  are  completely 

determined  by  the  given  pulse  p(t)  and  filter  basis  functions  ^(t). 

By  using  this  new  notation, 

Eqs.  (IV-36)  and  (IV-37)  become 

u(U  =  «jXa(t)  +  «2Xb(t) 

(IV-42a) 

x3(t)  =  -  »14>a(t)  -  o24>b(t) 

(IV-42b) 

As  before,  we  have  the  zero  bias  constraint  of  Eq.  (IV-16) 

rT 

*3(T>=J  p'(t)  x4(t)  dt  =  0 

{ IV —43) 

which  requires  that 

^  p'(t)  $a(t)  dt  +  ^  p'(t)  ^(t)  dt  =  0 

( IV -44) 

By  letting 

(T  p'(t)  *  (t)  dt 

pm  *-4? - 3 — 

/o  P'(t)  *b(t)  dt 

( 1V-45) 

Eq.  (IV -44)  then  requires  that 

«2  =  p(T) 

( IV— 46) 

and,  from  Eqs.  (IV-42),  we  have 

u(t)  =  Ixa(t)  +  p(T)  xb(t)]  <*i 

(IV-47a) 

xt(T)  =  -(4>a(t)  +  p(T)  *b(t)]  a4  . 

(IV-47b) 

The  minimum  cost  is  given  by 

-x42(T)/x5(T) 

( IV-48) 

where  now 

rT 

x4(t)  =  J  p"(t)  x,(t)  dt 

T 

=  -a1J  p"(t)  [*a(t) +  p(T)  *b(t)J  dt 

(IV-49a) 
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u2(t)  dt 


x5(T)  = 


=  <* \  j  lxa(t)  +  P(T)  xb(t)]2  dt  .  (IV-49b) 

The  minimum  cost  is  therefore 

if?  P"(t)  [*  (t)  +  p(T)  *.(t)]  dt}2 

pjr—  2  (IV-50) 

fQ  (Xa(t)  +  p(T)  xhlt)r  dt 


which  is  again  independent  of  the  unknown  multiplier  a Therefore,  we  can  set  =  1,  in  which 
case  the  optimum  solutions  are 

U(t)  =  xa(t)  +  p(T)  xb(t)  (IV-51a) 

Xi<t)  =  -*a(t)-p(T)  *b(t)  .  (IV-51b) 

More  importantly,  from  Eq.  (IV-46)  and  the  fact  that  a ^  =  1,  we  have  =  p(T).  Then,  from 
Eqs.  (IV-41), 

\4/2\5  =  1  (IV-52a) 

X3/2X5=p(T)  .  (IV-52b) 

Finally,  in  Eq.  (IV-35),  the  optimum  weights  are 
c  =  a  +  p(T)  b 

where  a,  b,  and  p(T)  are  defined  in  Eqs.  (IV-34a),  (IV-34b),  and  (IV-45),  respectively. 

In  effect,  what  we  have  done  is  to  project  the  matched  filter  onto  the  subspace  spanned  by 
the  set  of  functions  where  the  projection  has  been  done  with  respect  to  a  norm  involv¬ 

ing  the  mean-square  estimation  accuracy  of  the  filter  processor. 

From  this  example,  we  can  conclude  that  the  methodology  to  be  used  when  even  more  strin¬ 
gent  physical  realizability  constraints  are  to  be  imposed  is  relatively  straightforward. 

Example  3:  Resolution  constraints,  unrestricted  U 

In  this  case,  the  system  is  the  same  as  for  Example  1  with  the  additional  sidelobe  con¬ 
straints,  and  the  steps  in  the  synthesis  are  quite  similar.  For  this  reason,  the  details  are  pre¬ 
sented  in  the  Appendix.  Here,  it  was  not  possible  to  obtain  a  completely  analytical  solution  for 
the  optimum  control  because  of  the  large  number  of  constraints  which  had  to  be  satisfied.  How¬ 
ever,  we  were  able  to  reduce  the  variational  problem  to  a  nonlinear  programming  problem  in¬ 
volving  a  set  of  n  +  3  unknown  multipliers.  It  should  be  comparatively  easier  to  generate  the 
solution  to  the  optimization  problem  in  n  +  3  space  rather  than  to  deal  directly  with  the  general 
variational  problem  with  which  we  started. 

V.  CONCLUSIONS 

The  purpose  of  this  report  has  been  to  outline  the  methodology  to  be  used  in  formulating 
the  mismatched  filter  design  problem.  An  expression  for  the  range  accuracy  using  such  a  filter 
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has  been  derived  to  take  into  account  the  estimate  bias.  In  addition,  the  resolution  requirements 
of  the  radar  have  been  included  by  imposing  constraints  on  the  sidelobes  of  the  compressed 
pulse.  State-space  methods  were  then  used  to  develop  an  optimal  control  problem  whose  solu¬ 
tion  led  to  the  impulse  response  which,  in  turn,  would  yield  the  best  mean-squared  range  accu¬ 
racy  subject  to  preassigned  sidelobe  constraints.  Ignoring  the  sidelobes  and  applying  the  maxi¬ 
mum  principle  led  to  the  classical  matched  filter.  Control  techniques  have  been  used  before  to 
obtain  the  matched  filter  solution,  but  the  optimization  has  always  involved  a  rather  artificially 
defined  signal-to-noise  ratio.  Our  approach  takes  the  range  estimation  performance  into  ac¬ 
count  directly  and  exhibits  an  "optimum"  blend  of  the  communication  and  control  systems 
disciplines. 

Another  example  was  studied  in  which  the  sidelobes  were  again  ignored,  but  this  time  the 
set  of  admissible  filters  was  restricted  to  a  linear  combination  of  specified  basis  functions. 

The  optimum  weights  could  be  found  analytically  and  these  led  to  an  optimum  filter  within  a  sub¬ 
class  of  filters  not  necessarily  containing  the  matched  filter.  Such  a  realization  can  be  obtained 
using  RC  structures  which  are  relatively  simple  to  build. 

Applying  the  maximum  principle  to  the  design  problem  with  sidelobe  constraints  led  to  a 
nonlinear  programming  problem  involving  the  set  of  n  +  3  unknown  multipliers.  Implementing 
algorithms  for  solving  this  subsidiary  optimization  problem  should  be  easier  than  dealing  with 
the  original  variational  problem. 

In  the  most  general  problem  in  which  sidelobes  are  taken  into  account  and  in  which  the  class 
of  admissible  filters  is  restricted  by  physical  considerations,  numerical  methods  will  have  to 
be  used.  However,  computational  techniques  for  obtaining  the  optimal  solutions  are  well- 
established  in  the  control  field.  The  state-space  formulation  for  the  filter  design  therefore 
leads  to  a  technique  for  generating  the  optimum  mismatched  filter. 
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APPENDIX 


We  shall  now  derive  the  necessary  conditions  for  the  optimum  mismatched  filter  designed 

7 

for  sidelobe  reduction  as  well  as  for  good  range  accuracy.  In  this  case,  the  state  equations 
are  precisely  those  derived  at  the  end  of  Sec.  Ill  and  are  summarized  below  for  convenience. 


x^t)  =  -u(t) 

o 

II 

o 

X 

X2(t)  =  p(t)  X3(t) 

x2(0)  =  0 

x3(t)  =  p'(t)  X| ( t) 

x3(0)  =  0 

x4(t)  =  p"(t)  x3(t) 

x4(0)  =  0 

x5(t)  =  u2(t) 

Xg(0)  =  0 

Yj (t)  =  p(t  +  Tj )  x3(t) 

yj(0)  =0  j  =  1,  2,  .  .  .  , n 

(A-l) 

x3(T)  =  0 

y.2(T)  -  e.x2(T)4  0 

j  =  1,  2,  .  .  .  ,  n 

(A-2) 

Subject  to  the  preceding  conditions,  we  wish  to  find  the  control  function  u(t)  which  minimizes 
the  functional  —  x^(T)/xj.(T).  Notice  that  scaling  the  control  function  by  an  arbitrary  constant 
has  no  effect  on  the  terminal  constraints  or  the  cost  functional.  Therefore,  an  equivalent  prob¬ 
lem  is  to  minimize  —  x^(T)  subject  to  the  constraints  of  Eqs.  (A-l)  and  (A-2)  and  the  additional 
constraint 

x  (T)  =  1  . 

( A—  3 ) 

As  before,  we  define  the  functions 
9o[x(T)]  =  -x42(T) 
ejxfT)]  =  x3(t) 

02[x(T)]  =  xg(T)  -  1 

^(x<T),y<T)]  =  y2(T)  -  e.x2(T)  j  =  1,  2 . n  .  (A-4) 

In  this  case,  the  Hamiltonian  is 

H[x(t),  y(t),  u(t),  X_(t),  p(t)]  =  — A1(t)  u(t)  +  X2(t)  p(t)  x3(t)  +  X3<t)  p(t)  Xj(t) 

+  X4(t)  p"(t)  Xl(t)  +  X5(t)  u2(t) 
m 

+  Z  p(t  +  t  )  y.(t)  .  (A-5) 

j  =  l 
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The  costate  variables  satisfy  the  equations 


yt) 


OH 

8xk(t) 


which,  in  this  case,  are  given  by 


( A- 6) 


Ad(t)  =  -p(t)  A2(t)  -  p1  (t)  A3(t)  -  p"(t)  A 4( t )  -  P<t  +  Tj)  M-j(t) 

j=1 


Ak(t)  =  0 


k  =  2,  .  .  .  ,  5 


4j(t)  =0  j  =  1,  2,  .  .  .  ,  n 

The  terminal  values  of  the  costate  variables  are  obtained  from 


1  36  |x(T)l  "  8*.[x(T),  ^ITl] 

VTI  ■  2  “i  -talrT  *  2  -i  -1 - 

J  J 


i=0 

n 


i=l 


dx^T) 


_  dip.[x{T),  Z(T)] 

lVT)  =  Z  »,  -'-avIT) 

i=l  3 


and,  therefore, 

Ad(T)  =  0 


A4(T)  =-2aox4(T) 


X2(T)  =-2x2(T)  £ 
i=l 


Mj(T)  = 


j  =  1,  2,  ....  n 


A3(T)  = 


In  addition,  the  constant  multipliers  are  zero  if 


y.2(T)  -  e.x^(T)  <  0 


(A-7a) 


(A-7b) 

(A-7c) 


(A-8a) 


(A-8b) 


(A-9) 


and  v.  ^  0  otherwise.7  Therefore,  we  assume  that  the  summations  on  j  include  only  those  terms 
corresponding  to  sidelobe  constraints  which  are  actually  on  the  boundary.  Now  we  combine 
Eqs.  (A-7)  and  (A-9)  to  get 


n 


\2(t)  =  -2x2(T)  £  S^i“X2 

A5(t)  =  a2  =  A 

i=l 

A3(t)  =  a1  =A3 

fi.(t)  =  2y.(T) 

A4(t)  =  -2a?ox4(T)  =  A4 
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and  the  equation  for  X^(t)  is  simply 


X1(t)  = -X2p(t)  -  X3p'(t)  -  X4p"(t)  -  2  (Jup(t  +  Tj) 

j=l 

with  X^(T)  =  0.  This  equation  is  easily  integrated  to  give 

r»T 


(A-ll) 


j=1 

As  before,  we  assume  that 

p(0)  =  p(T)  =  0 
p '  ( 0 )  =  p'(T)  =  0 

so  that 


Xi(t)  =X2  J  p(cr)  d(T  +  X3[p(T)  -  p(t)]  +  X4[p'(T)-p’(t)] 

n  rp 

Z  ^  \  p(tr  +  d<T  • 

J  J 


(A- 12) 


(A-13a) 
( A-l 3b) 


Xj(t)  =  —  X4p(t)  -  X3p(t)  +  X2  ^  p(a)  dir  +  Z  §  p(ff  +  t  )  da  (A-14) 

‘  j  =  l  1 

The  maximum  principle  states  that  the  optimum  control  must  maximize  the  Hamiltonian. 
Therefore,  it  is  necessary  that  the  function 


h(u)  =-X1(t)  u(t)  +  X5u^(t) 


(A-15) 


be  maximized  for  the  optimum  u  e  U.  We  shall  assume  that  U  is  arbitrary  and  that  a q  <  0, 

<  0.  (This  must  be  the  case  from  arguments  given  in  Sec.  IV.)  Therefore,  the  optimum  con¬ 
trol  is  of  the  form  u(t)  =  X.(t)/2X_.  Using  Eq.  (A-14),  we  conclude  that 

1  D 


S(t)  =  2X^  P'(t)  -  2j^  P(t)  +  ^  P(a)  dtr  +  Z  P(<T  +  Tj)  d(T  •  (A 


16) 


a  and  z(t): 


to  use  vector  notation  to 

express 

~V2X5 

Zjtt)  = 

p'(t) 

-X3/2X5 

z2(t)  = 

p(t) 

+X2/2X5 

z3(t)  = 

r  p(.,  . 

jt 

PJ/2X5 

zj  +  3(t)  = 

f  * 

da 


J  da  j  =  1,  2, 


(A- 17) 
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The  optimum  control  may  then  be  written  as 


u(t)  =  Qftz(t) 


(A-18) 


where  the  superscript  t  denotes  the  vector  transpose.  We  define  another  n  +  3  vector  u;(t)  by 
setting 


—(t)  ~  ~  z(<r)  da 


Then,  from  Eqs.  (A-l),  the  optimum  filter  is 
x1(t)  =  t) 


(A- 19) 


(A-20) 


It  is  now  possible  to  solve  for  the  terminal  values  of  the  remaining  state-variables  by  substitut¬ 
ing  in  Eqs.  (A-l) 

"T  .  x  nT 


nl  nl 

X2(T)  =  \  p(t)  x1dt  =  o1  \  p(t)  0)(t)  dt 
Jo  Jo 

rT  +  rT 

x3(T)  =  j  p'(t)  x1(t)  dt  =  J  p’(t)  a>(t)  dt 
T  r*T 

x4(T)=J  p"(t)  xddt  =  a}  j  p"(t)  w(t)  dt 

x5(T)  =  ^  u2(t)  dt  =  dtJ  £ 

rT  +  rT 

y.(T)  =  \  p(t  +  r.)  x.(t)  dt  =  a/  t  p(t  +  r.)  w(t)  dt 
j  J  1  J 


'’T  -  .  rT 

^o  J  ^o 

It  is  convenient  to  introduce  some  additional  notation: 

n  T 

Jo 

"T 


(A-21) 


Zpw  =  §Q  P(t)  ^(t)  dt 

Ip'w=  Jo  P’(t)^(t)dt 

vw-r 


pM(t)  w(t)  dt 


“O 

"T 


y .  =  C  p(t  +  t.)  u>(t)  dt 
-j  J 


(A-22) 


where  all  the  above  vectors  have  n  +  3  components.  In  addition,  we  define  the  (n  +  3)  x  (n  +  3) 
matrix 


=  ^  z(t)  2} 


(t)  dt 


(A-23) 
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so  that 


*2(T) 

t 

=  a  y 

—  -pw 

x3(T) 

t 

=  a  y  , 

-  ip'w 

x4(T) 

=  a  S' 

—  -p’  rw 

x5(T) 

=  aS* 

y3(T) 

ii 

& 

It  is  important  to  note  that  tne  vectors  y  and  matrix  r  are  completely  known  quantities.  Only 
the  vector  a  needs  to  be  determined,  and  it  must  satisfy  the  constraints  of  Eqs.  (A-2)  and  (A- 3), 
namely, 

x3(T)  =  0 

y.2(T)-ejx22(T)4  0 

x5(T)  =  1  j  =1,2 . n  .  (A-25) 

These  become 


o  y  ,  =0 

-  -pw 

t  t  t  .  n 

a  y  .y .  ct  —  «: .  a  y  a  <  0 

—  -J-3  —  J-  ipw-^ 


a  T  a  =  1 


(A-26) 


These  equations  represent  n  +  2  constraints  on  the  n  +  3  variables 
degree  of  freedom  is  used  1o  minimize  the  cost  functional  —  x^(T), 


off,  .  .  .  ,  an+y  The  remaining 
which  in  this  case  is 


Now,  we  let 


t  I 

—  ot  y  ,,  y  ,,  a 
—  -Lp' '  w— p  1  w  — 


r .  =  y.y  }  —  e  .y  y  * 

J  -J-J  jlpwipw 


( A-27 ) 


r0  ■  V'w-p"w) 

and  the  original  variational  problem  for  the  optimum  mismatched  filter  reduces  to  the  following 
nonlinear  programming  problem.® 

Find  the  n  +  3  vector  a  which  maximizes  the  function 


f(  a)  =  a*r  a 
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subject  to  the  constraints 


rAa)  =  a'y  =  0 
1  —  -pw 

r  a)  =  g^Ta  —  1  =  0 

j  =  1.  2 . n  .  (A-28) 

J  J 

It  may  be  that  this  nonlinear  programming  problem  is  easier  to  solve  than  the  original  variational 
problem. 
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